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Abstract 

It is shown that the time evolution of the squeezed and displaced state 
may be obtained by solving the Heisenberg equation of motion of an appro- 
priate operator and finding the eigenstates of the time evolved operator. The 
connection between symplectic transformations and squeezing is explored. 
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1 Heisenberg equation of motion 

It is well known that the squeezed and displaced states of the harmonic oscillator 
are eigenstates of a linear combination of the creation and annihilation operators 
of the oscillator with non-zero eigenvalue (Yuen 1975, Henry and Glotzer 1988). 
The properties of squeezed states have been discussed extensively,see for example: 
Ekert and Knight 1989, Gong and Aravind 1990, Gersch 1992. The squeezed state 
is a limiting case of the squeezed and displaced state corresponding to eigenvalue 
zero. Since the Heisenberg equations of motion for the position and momentum 
operators of a simple harmonic oscillator can be solved, it is possible to obtain 
the time evolution of the squeezed and displaced state by solving the Heisenberg 
equation of motion for a linear combination of position and momentum operators 
and then finding the eigenstates of the time evolved operator and establishing a 
relation between the two states. We outline such a procedure and also show that 
there is a connection between squeezing and symplectic transformations of position 
and momenta. For simplicity of notation, Planck's constant h, the mass parameter 
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m and the frequency of the oscillator oj are all set equal to 1. In the first part of 
this paper we establish a relation between the eigenstates of a Heisenberg operator 
and the solutions of the time dependent Schrodinger equation. We then use this 
relation to obtain the time evolution of the squeezed and displaced state of the 
harmonic oscillator. In the final part of the paper we consider the connection 
between squeezing and symplectic transformations. 

Consider the eigenvalue equation for an operator A given by 

A(t) = \j(t) (1) 

The Heisenberg equation of motion for an operator A which has no explicit time 
dependence is given by 

f)A 

'^ = IAH] (2) 

where if is a time independent Hamiltonian. Equation (2) has the formal solu- 
tion 

A(t) = e iHt A(0) e~ iHt . (3) 
The eigenvalue equation at time t=0 given by 

A(0)^(0) = A^O) (4) 

may be rewritten using the solution given by equation (3) in the form 

A(t) [e im $,(0)] = A, [e iHt $,(0)] (5) 

which shows that the eigenvalues Aj of A(0) are also eigenvalues of A(t) and the 
eigenstates of A(t) are given by 

~ e im $,(0). (6) 

In the above equation the ~ symbol is used because a comparison of equations (1) 
and (5) shows that it is possible for the two normalised states to differ by a phase 
factor which depends only on time but not any other variables that if is a function 
of, as will be demonstrated explicitly in the next paragraph. Since the eigenvalues 
are time independent we can make the identification that Aj(t) in equation (1) 
has no dependence on t, i.e. Aj = Aj(t). It may not always be convenient to 
solve equation (4) and then use equation (6) to find the eigenstates at time t as 
$j (0) may be a complicated linear superposition of eigenstates of the Hamiltonian 
and the evaluation of the right hand side of equation (6) in closed form may be 
difficult. Under certain circumstances it may be possible to find A(t) directly from 
equation (3) and solve the eigenvalue equation for A(t) given by equation (1) to 
find the eigenstates and eigenvalues of A(t). The equation satisfied by $j(t) may 
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be found by differentiating equation (f) with respect to time and using equation 
(2) to give 

[A -A,} [H^(t) + = 0. (7) 

This equation is clearly satisfied if 

[H + i^] = N (t) *,.(*) (8) 

in which Hj(t) is an arbitrary function of time. The solutions of equation (8) may 
be given in the form 

^j(t) = e-'finW'* Fj(t) (9) 

where 

[H + Fjit) = 0. (10) 

It is now possible to identify that 

= Fji-t) (11) 

is the Schrodinger state evolving under the influence of H from the state <&j (0) = 
Fj(0) = ^j(O). Thus we can relate the Schrodinger state ^j(t) to the solutions of 
equations (1) and (3) in the form 

^.(t) = e- <Ht *j(0) = e im $j(-t) (12) 

where ^ 

5(t) = [ N (t') dt'. (13) 



We have shown that apart from a time dependent phase factor, the solution ^ j{t) 
of the Schrodinger equation for H may beconstructed from the solution <S>j(—t) 
of the eigenvalue equation for the time dependent Heisenberg operator A(—t) by 
propagating backwards in time. To find the phase factor S(t) we can expand both 
*&j(t) and &j(t) in terms of the eigenstates of H with eigenvalues E n which are 
denoted by |n), n — 0, 1, 2, .... Let 

*,-(t) = $>„(t)|n>, 

n 
n 

It is clear from equations (12) and (13) that 

a n (t) = b n (-t) e iS W, 

On(0) = 6 n (0) (15) 
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In terms of the eigenvalues E n of H 

a n (t) = a n (0) e~ iEnt . (16) 

Hence 

e i 5(t) = -iE n t b n(0) / 17 x 

bn(-t)' K ' 

The simplest way to find 5(t) is to choose n = so that 

e , S( „ = e -,« w±m (18) 

(0|*i(-*)> 

Thus the Schrodinger state ^j(t) may be obtained from the solution $.,•(£) of 
equations (1) and (3) using the relation 

= e - i£ot ^ $ (19) 

2 Time evolution of squeezed states of simple 
harmonic oscillator 

We now illustrate the procedure by studying the time evolution of the squeezed and 
displaced state of the simple harmonic oscillator. Consider the harmonic oscillator 
governed by the Hamiltonian 

H = -{P 2 + X 2 ). (20) 

The commutation relation [X, P] = i maybe used to establish the Heisenberg 
equations of motion for the operators X(t) and P(t) which may then be solved to 
give 

X(t) = X(0) cost + P(0) sint, 

P(t) = P(0) cost - X(0) sini (21) 

Let 

A(t) = iaP(t) + px(t). (22) 

Then using the operator solutions given by equation (21) and adopting the notation 
that X(0) = X, P(0) = P, a(0) = a, /?(0) = /?, A(t) may be given in the form 

A(t) = ia(t)P + p(t)X, (23) 
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where 

a(t) = a cost — i/3sint, 

fj{t) = p C ost - zasint. (24) 

Using P = i-J^£ it is easy to verify that the eigenstates of A(t) with eigenvalue A 
are given by 

$(t) ~exp^ P{t) <v A 



( x ~m) 



2a(i)V" /3(t)' (25) 

It is clear that $(0) represents a squeezed and displaced state with squeeze pa- 
rameter S and displacement parameter D given by 

5(0) = £ , D(0) = £ . (26) 

a p 

The time evolution of the squeezed and displaced state may now be obtained using 
equations (19), (25) and (24). It is easily seen that the squeeze and displacement 
parameters of the time evolved state are given by 

/3(—t) (/3 cost + zasint) 
a(—t) (a cost + zpsint) 

D(t) = = , -j , , y (27) 

P{~ t) (pcost + zasmt) 

These equations may be expressed in terms of the squeeze and displacement pa- 
rameters at time t = in the form 

(5(0) cost + zsint) 
(cos t + iS(0) sin t) ' 

(o(O)cost + zsmtj 

The squeezed and displaced states are displaced gaussians with complex exponents. 
The normalisation factor associated with these states and the overlap integrals 
appearing in equation (19) depend upon t and may be readily evaluated. The 
normalised squeezed and displaced state may finally be given in the form 

V sd (x,t) = N(t) exp[- ls(t)(X 2 - 2D(t)X + D(t)D(0) cost)] , 

N(t) = (^®)* (cost + iS(0) sin t)~* . (29) 
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The squeezed vacuum state corresponds to choosing D(0) = and coherent states 
correspond to choosing 5(0) = 1 in the general expression given above. It may be 
shown that the probability distribution associated with ^ s d(t) is given by 



^ sd = (^exp[- 7 (X ■ D(O)cosi) 2 ], 



1 (cos 2 t + S 2 (0) sin 2 t) ' { } 

This probability distribution corresponds to a wavepacket whose centre oscillates 
with classical frequency 1 between the classical turning points ±D(0) and the 
shape of the wavepacket undergoes a time dependent squeeze. The simple feature 
that is a gaussian with time dependent squeeze and displacement parame- 
ters suggests that there may be a connection between squeezing and canonical 
transformations in classical mechanics. Towards this end we consider a linear 
transformation from a set of canonical coordinates X\ and P\ satisfying the com- 
mutation relation [Xl,Pi] = % to a new set of coordinates X 2 and P 2 through a 
mapping of the form 

X 2 = aX 1 + bP-y , 

P 2 = cX x + dP 1 . (31) 

It is clear that if 

ad -be = 1 , (32) 

then [X 2 ,P2] = i and the new coordinates also satisfy the same commutation 
relation as the old coordinates. Transformations of the form of equation (31) with 
unit determinant are symplectic transformations (Goldstein 1980). It is well known 
that apart from the unit matrix there are only three other matrices of dimension 
2 which have unit determinant (Schiff 1968) and they may be written in terms of 
the Pauli spin matrices 

°* = ( i o ) ' ff y = ( +i ~o ) ' ° z = ( o -l ) (33) 

in the form 

_ e<Tx _ f cosh^ sinh^ \ ivCy f cos// sinz/ \ ( e p 



e _ Q°-"y — e - 

y sinh 9 cosh 9 J ' y — sin v cos v J ' \ e~ p 

(34) 

Any matrix in two dimensions with unit determinant may thus be constructed 
from the exponentials of the matrices l,a x ,a y and a z . If we now consider the 
eigenvalue equation 

(iaPi + pX^ = (35) 
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then 

= exp(-^ 1 X 1 2 ) , S 1 = ^ (36) 

is a squeezed state with squeeze parameter Si. In terms of X 2 and P2 equation 
(35) may be written in the form 



with the solution 



i(aa + i(3b)P 2 + (/3d - iac)X 2 V = (37) 

* - • s * ^ ^Tli • (38 > 

which is also a squeezed state with squeeze parameter S 2 . Thus we can see that 
when two sets of coordinates are connected by a symplectic transformation then 
a squeezed state with respect to the old coordinate transforms into a squeezed 
state with respect to the new coordinate with a simple transformation between 
the squeeze parameters given by 

S 2 = S 4^{. (39) 

Comparison of equations (21) and (31) shows that 

a = cost, b = sini, c = — sini, d = cost. (40) 

Since (ad — be) = 1 it is clear that the time evolution of the position and momen- 
tum operators given by equation (21) is a symplectic transformation. The time 
evolution of the squeeze parameter given by equation (28) is of the form given 
by equation (39). Comparison with equation (34) shows that time evolution is a 
symplectic transformation of the type induced by a y . In this paper we have out- 
lined a novel procedure for studying the time evolution of the Schrodinger states. 
We have shown that the time evolution of the squeezed and displaced state may 
be studied using this method. We have shown that under a symplectic transfor- 
mation a squeezed state transforms into a new squeezed state with a transformed 
squeeze parameter. For a 1-dimensional system there are three sets of symplectic 
transformations possible apart from the identity transformation. We have shown 
that the time evolution of the squeezed state may be related to one of the three 
non-trivial symplectic transformations. 
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